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Abstract: In this paper we study possible observational consequences of the 
bouncing cosmology. We consider a model where a phase of inflation is preceded 
^ \ by a cosmic bounce. While we consider in this paper only that the bounce is due 

to loop quantum gravity, most of the results presented here can be applied for 
different bouncing cosmologies. We concentrate on the scenario where the scalar 
field, as the result of contraction of the universe, is driven from the bottom of 
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O \ the potential well. The field is amplified, and finally the phase of the standard 

slow-roll inflation is realized. Such an evolution modifies the standard inflation- 
ary spectrum of perturbations by the additional oscillations and damping on the 

k> \ large scales. We extract the parameters of the model from the observations of 

the cosmic microwave background radiation. In particular, the value of inflaton 
mass is equal to m = (2.6 ± 0.6) ■ 10^^ GeV. In our considerations we base on the 
seven years of observations made by the WMAP satellite. We propose the new 
observational consistency check for the phase of slow-roll inflation. We investi- 
gate the conditions which have to be fulflUed to make the observations of the Big 
Bounce effects possible. We translate them to the requirements on the parame- 
ters of the model and then put the observational constraints on the model. Based 
on assumption usually made in loop quantum cosmology, the Barbero-Immirzi 
parameter was shown to be constrained by 7 < 1100 from the cosmological ob- 
servations. We have compared the Big Bounce model with the standard Big 
Bang scenario and showed that the present observational data is not informative 
enough to distinguish these models. 
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1. Introduction 

The observations of the cosmic microwave background (CMB) radiation indicate 
that the power spectrum of primordial scalar perturbations is in the broad range 
nearly scale-invariant. Therefore, the spectrum can be written in the power-law 

form 

k ^ "=-' 



V.ik) = Ai^-j , (1.1) 

where the spectral index Ug is close to unity. Here, Ag is an amplitude of the scalar 
perturbations and ko is the so-called pivot number. The case Ug = I corresponds 
to the scale-invariant Harrison-Zeldovich spectrum. The observed spectrum is 
almost of this type, namely it is little red-shifted (rig < 1). In particular, the 
seven years of observations made by the WMAP satellite |]l[ indicate that Ug = 
0.963 ±0.012 (68% CL). 
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The spectrum in the form discussed above can be explained by the phase of 
cosmic inflation (see e.g. [^). This phase can be driven by the self- interacting 
scalar field, the so-called inflaton field. In the most conservative approach the 
inflation can be driven by a single massive scalar field. This case will be consid- 
ered in this paper. In this model, a nearly scale-invariant spectrum spectrum is 
generated during the slow- roll phase. In a more general case the inflation can be 
driven with the different potentials. However, the potentials other than massive, 
lead to the non-Gaussian structure of the cosmic primordial perturbations what 
can be constrained by the CMB observations [Q. Since the CMB anisotropics 
do not indicate any non-Gaussian signatures, the massive potential is somehow 
privileged. However, with the present sensitivity on these kind of effects, some 
of the other potentials are still allowed. Other models, as multi-field inflation 
are also possible to be realized. In this paper, we consider the simplest possible 
realization of the inflation which is given by the single massive scalar field. 

The weak point of the slow-roll inflationary scenario is that it requires some 
special initial conditions. Namely, the field has to start its evolution not from 
the bottom of potential well but from the position which is far from its center. 
In the classical model there is no mechanism to drive this field up the potential 
well. However, it has been recently pointed out [|| that the phase of a quantum 
bounce can drive the inflaton field up the potential well and set the proper initial 
conditions for the slow-roll phase. The studies were performed within the loop 
quantum cosmology (LQC) p, ^, however the mechanism is generic for all models 
with the bouncing phase. In the framework of LQC this issue has been studied 
recently in 1,0, 1, 1 . 

In the framework of LQC the classical dynamics of the universe is significantly 
modified when the energy density approaches the Planck energy density. These 
effects of the quantum gravitational modification can be introduced as correc- 
tions to the classical equations of motion. In particular, the modified Friedmann 
equation takes the form 

where the critical energy density is defined as follows 

where ppi := m-pj and mp\ ~ 1.22 ■ 10^^ GeV is the Planck mass. The above 
expression ( p..3|) is crucial, because it relates the Barbero-Immirzi parameter 7 
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with the parameter pc, which can be constrained observationally. This is a good 
example of expressing the phenomenological parameter, as pc, in terms of pa- 
rameters of the underlying theory. In this case the fundamental theory is loop 
quantum gravity (LQG) with a free parameter 7. 

The value of Barbero-Immirzi parameter is usually fixed from the considera- 
tions of the black hole entropy [|10|- In particular, the value 7 = 0.239 ||lT| was 
derived what leads to pc = 0.82 m-pj. The expression (|1.3|) is however not free 



from ambiguities. In particular, it bases on assumption that the area of the loop 
in LQC is equal to a gap of the area operator within LQG. This is not necessary 
true and therefore findings based on expression ( |1.3| ) must be carried with a due 
care. In general, pc can be treated as a free phenomenological parameter (see e.g. 



12, 13, 14 



2. Inflation in LQC 



The global dynamics of the considered model was studied in |T5|. It was shown 



there that the model possesses generic inflationary attractors. However, the par- 
ticular evolutionary paths can differ. Here we restrict our considerations to the 
one particular scenario where the energy density is dominated by the kinetic term 
at the bounce. This condition guarantee that the quantum back-reaction effects 
can be neglected and equation ( |1.2| ) can be applied. This issue will be discussed 
in more details at the end of this section. 

In the considered evolutionary scenario, we initiate evolution from the con- 
tracting phase. The scalar field is initially placed at the bottom of the potential 
well, what seems to be the realistic and conservative assumption. The initial 
energy density is therefore contained in the kinetic part only. Because we begin 
the evolution at the low energy scales, the corresponding time derivative of the 
field is also small. Another possible choice of initial conditions is given e.g. by 
saturating the Heisenberg uncertainty relation, then vr^(to)0(^o) = ^/2 at some 
initial time to- 

The field begins its evolution from some tiny quantum fluctuations and starts 
to oscillate at the bottom of the potential with the time scale of oscillations 
proportional to m^^. The dynamics of the field is governed by the unmodified, 
standard equation 

+ 3i70 + mV = O. (2.1) 

During the contraction, the Hubble factor is negative, and therefore the second 



term in equation (|2.1|) acts as anti-friction. Moreover, while the universe contracts 



3 



the absolute value of H increases. Therefore as an effect of growing anti-friction, 
the oscillations are amplified. While approaching the bounce, the field becomes 
rapidly displaced from the equilibrium state. After that, the standard slow- roll 
inflation starts. We present the described evolution of the scalar field in Fig. ^ 



The described evolution cor- 
responds to the shark fin sce- 
nario discussed in [|]. During 
the evolution, the scale factor 
decreases in the pre-bounce stage 
Since in this stage the field be- 
haves effectively as a dust mat- 
ter, the scale factor a oc |tp'^. 
The scale factor reaches its min- 
imal value and then expands 
almost exponentially, during the 
slow-roll inflation. While the 
inflation ends, the scale factor 
increases as a oc \t\^/'^. The 



4' 

3- 
2- 
1 

0- 
-1 




-40 



-20 







20 



40 



mt 



Figure 1: The shark fin type evolution of the scalar 
field for m = 10~^mpi. Here 0max ~ 3?7ipi, what 
leads to the total number of e-foldings A^ ~ 56. 



numerically computed evolution of the scale factor has been shown in Fig. |[ 

We finish this section with 
some remarks on the quantum 
back-reaction effects. These ef- 
fects become significant when 
the evolution of the higher mo- 
ments influences behavior of the 
mean values of the quantum 
operators. In the free field case, 
the evolution of the mean val- 
ues decouples from the evolu- 
tion of the higher moments (for 
detailed discussion of this is- 
sue we refer to [|16|). The back- 
reaction effects appear however 
in the presence of the poten- 
tial of the scalar field and can significantly modify the effective dynamic. In 
particular, equation (|1.2|) holds only when the quantum back-reaction effects can 




Figure 2: Evolution of the scale factor for m 
IQ-'^mpi. 
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be neglected. Otherwise this equation can be generahzed to |T7| 



1 da\' 
a dt J 



6TT 



3mpj 



Pc 



± 






(2.2) 



and the dynamics can be much more comphcated that this discussed previously. 
Here 77 parameterizes strength of the quantum back-reaction effects. In the con- 
sidered massive field case the energy density and pressure are respectively 



p = 1 and F = . 

'^ 2 2 2 2 



(2.3) 



bounce 
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When P = p, the quantum back-reaction effects disappear, equation ( p.2| ) simplify 
to (|1.2|) . This corresponds to the free field case. Therefore, while the energy 
density is dominated by the kinetic part, the quantum back- reaction effects can 
be neglected. In Fig. ^ we show an exemplary evolution of the energy density 
in the considered model. We also present contributions from the kinetic and 
potential parts. 

The kinetic term ^ dom- 
inates the potential part ^^^^ 
in the broad region around the 
bounce. This is the shadowed 
region from ti to t^ in Fig. ^ 
Therefore approximation based 
on (|1.2|) holds. For the densi- 
ties p <^ Pc the kinetic part 
can be dominated by the po- 
tential part. However, at these 
densities the quantum effects 
become unimportant and dy- 
namics can be approximated 




slow-roll 
inflation 



Figure 3: The solid (black) curve represents the 
evolution of the total energy density of the scalar 
field. The dashed (blue) curve represents the con- 

by the classical equations. There-*"^^^^"^ ^^°™ *^« ^^"^*^': P^^*' ^he dotted (red) 



fore condition ^ ^ 



m^tf , ,1 curve represents the contribution from the potential 
part. In the filled region the energy density is dom- 
inated by the kinetic part. Here we have assumed 
10-4. 



m 



^mpi. 



be fulfilled only in the vicinity 

of the bounce. If it is not, the 

quantum back reactions must 

be taken into account. At time ^2? the bounce takes place and the energy density 

reaches its maximal value pc- Later, the energy density decreases and holds at 

approximately constant value. This is the sign that the phase of inflation starts. 
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Precisely, the inflation starts at time t^ when the field turns round. At this point, 
the kinetic term falls to zero since the field stops for a moment. Thereafter, the 
field makes a slow-roll from the top of the potential well. This part of evolution 
is almost purely classical and the quantum corrections can be neglected. 

3. Observational hints on the slow- roll inflation 

Before we proceed to investigate the possible effects due to the bounce, we will 
firstly discuss the present observational hints regarding the infiation. It is crucial 
since probing the infiationary phase is more observationally available and give 
the chance to fix some parameters of the model. We will discuss here what we 
can already say about the slow-roll infiation model in light of the latest WMAP 
observations. We also propose the consistency check on the slow-roll infiationary 
scenario. 

The seven years of observations made by the WMAP satellite give the fol- 
lowing values of the amplitude and spectral index of the scalar perturbations 

i 

A = 2.441+°;°^^ 10-^ (3.1) 

ns = 0.963 ±0.012, (3.2) 

at the pivot scale ko = 0.002 Mpc~^ (see equation ( [L.l| )). The prediction from 
the slow-roll is the spectrum of scalar primordial perturbations in the form 



nk) = - —] 777 , (3.3) 




as well as the spectrum of the tensor perturbations (gravitational waves) in the 

form 

16 / if \ V A; ^ "' 



Expressions for the scalar and tensor spectral indices are respectively 

ns = l + 2r]- 6e, (3.5) 

and 

rit = -2e, (3.6) 
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where e, r/ ^ 1 are called slow-roll parameters. The rj and e are defined in the 
following way 

ml, fV"\ _ml,l 



so for the massive infiation field, rj = e. Based on this and equation (|3.5| ) we find 
e = (1 - n,)/4. With use of Q, this gives us e = 0.010 ± 0.003. Therefore, the 
slow-roll condition e ^ 1 is indeed fulfilled. Moreover, based on (|3.2|) , the tensor 
spectral index is predicted to be 

n — 1 
rit = — = -0.019 ± 0.006. (3.9) 

Above we have related observations of the spectral index with the expression 
predicted from the slow-roll infiation. It was straightforward since both, predicted 
spectrum ( p.3|) and spectrum used in fitting (|1.1|) had the same power-law form. 
However relating the fitted parameter A^ with S requires additional discussion. 
At the pivot scale Vs{k = k^) = Ag. Moreover, we know that the infiationary 
spectrum at a given mode k is formed when this mode crosses the horizon, namely 
when k ~ aH. Afterward, the spectrum holds the form fixed at the horizon. 
Therefore at the given scale k, the spectrum is Vs{k = aH) = S. Since S decreases 
with time, the observed spectrum has the falling tendency (governed by the power 
law dependence). Based on this, one can relate Vs{k = k^) = Vs{k = aH), what 
gives As = S. Therefore, As gives us the value of S at the point when the mode, 
which is at present equal to ko, had crossed the horizon during the infiation. This 
observation will be crucial for the later considerations. 

In order to quantify the contribution from the tensor modes it is convenient 
to consider the ratio 

r ^ 'If^^l'l = I = 16e = 4(1 - n^) = 0.15 ± 0.05, (3.10) 

Vs{k = ko) S 



where in the last equality we have used the WMAP results (|3.2|). This result is 



consistent with the present constraints on the contribution from the tensor modes 

r < 2.1 at 95% CL (WMAP-7, 0), (3.11) 

r < 0.73 at 95% CL (BICEP, pi). (3.12) 
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Moreover, the value predicted in ( |3.10|) is above the observational threshold on 
detection of the PLANCK satellite [|^ . Therefore, if the predictions of the slow- 
roll inflation are correct, the tensor modes should be observed by the PLANCK 
mission. 

The next parameter that can be computed is the value of the scalar field. 
Combining (pl7|), (|3.8|) and (p.5|) we get 

mpi 



''obs 



y^ 



2.9±0.5mpi. (3.13) 



n. 



At this value of field, the observed structures were created ^. Therefore it can 
be treated as an lower limit on the maximal displacement of the scalar field. 
The maximal value of the scalar field is unbounded within the classical theory. 
However LQC puts the constrain on its value since the energy density is bounded 
by pc- Based on this one find that 

101 < ^^- (3.14) 

m 

The found value of the inflaton field (|3.13|) can be translated into the corre- 
sponding e-folding number 

AT^bs ^ 27r% = -^ = 54 ± 18. (3.15) 

mpj 1 — ^s 

This is also not the total e-folding number for inflation, but only the lower limit 
on its value. The total e-folding number is, in LQC, constrained by 

N<^^„ (3.16) 

what bases on (p.l4|) . Finally one can also derive the mass of inflaton field. 
Namely 

m ~ mpi-A/37rAs(l - ris) 

= (1.4 ±0.5) -lO^Wi 

= (2.6±0.6)-10^3GeV. (3.17) 

Therefore, one of the parameters of the model is fixed. The remaining parameter 
Pc is however harder to determinate. We will discuss the present observational 



^Precisely it is the value of 4> at which mode which is at present equal to fco, had crossed 
the horizon during the inflation. 



constraint on pc in Sec. |. We stress that it was possible to determinate the 
value of the parameter m basing only on the observational effects of inflation. It 
was not necessary to introduce any LQC effects here, because they were negligible 
during the phase of inflation. In other words, the slow-roll inflation is the classical 
(and observationally available) part of the considered shark fin scenario. 

The usual consistency check of the inflationary models bases on expressing 
of the tensor-to-scalar r in terms of others (measured) parameters of the model. 
In case of the slow-roll inflation, we have considered it in equation ( p. 101 ) and we 



have shown that the derived value of r places within the observational bound. 
Here, we propose an additional consistency check for inflation. This new consis- 
tency relation requires however information about the duration of the reheating 
phase. Alternatively, the method can be used to put a constraint on the phase 
of reheating after inflation. 

The consistency check base on the fact that the modes created at the par- 
ticular point of inflation where = 0obs, correspond to the present pivot scale 
at which the amplitude of perturbations was computed. These particular modes 
have the size of horizon when created from the quantum fluctuations. Therefore, 
at this particular point 

^H = 2vr^ ^ 1 = ^ = 1.2 • 10^/pi. (3.18) 

In turn, the present pivot scale is equal to 

27r 
Ao = — = 3.14 ■ 10^ Mpc. (3.19) 

Based on this, one can find the total increase of the scale factor from the point 
at which = 0obs, fiU now. We obtain the value 

Atot := ^ = ^ = 5 ■ 10^1 (3.20) 

There is also another way to compute this quantity. Namely, starting from 
= 0obs, the length Ah grows thereafter till the end of inflation, across the 
reheating, radiation domination phase, matter domination phase until now. It 
is hard to precisely determinate the increase of the scale factor at this whole 
evolution. In particular, because we do not know the duration of reheating phase 
and the duration till the end of inflation was determined with the significant 
uncertainty. Therefore we can perform only a raw approximation of the total 
increase of the scale factor. 
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Let us collect the particular contributions starting from the present and going 
backward: 

• Matter era. The period from recombination till now. A^at •= 1 + ^dcc — 10^ 

• Radiation era. The period from the end of reheating till the recombination. 

A , ._ Tout _ 10^* GeV ^ c . i r|23 
""^ '~ 'Tk^ ~ 0.2 eV - -• -LU 



• Reheating (see e.g. |20, 21|). The period when the particles are created from 
the decaying inflaton field and the universe thermalizes. The corresponding 
quantity Areh is model dependent and should be fixed for the particular 
scenario. For instance, for the instantaneous reheating Areh ~ 1- However, 
for the considered chaotic inflation the duration of reheating can be longer. 
We leave detailed considerations to this issue for the further studies, and 
now let the value of Areh as a free parameter. 

• Inflation. The increase of the scale factor is equal to Ai^f := 6^°*== ~ 3 ■ 10^^ 
Based on this, the total increase of the scale factor is equal to 

Atot = AinfArehAradAmat- (3.21) 



The left side in equation ( |3.21 ) is determined from ( p. 20 ). Based on the above 



relation one can e.g. try to determine duration of reheating. Namely, we have 



At,ot 5 ■ 10^5 



^reh 



Ai„fAradA„,at 1.5 ■ 1051 



3 ■ 10^ (3.22) 



On the other hand, in order to use ( p.21| ) to verify the model of inflation, the 



duration of reheating must be known from the theory. Then one can define the 
quantity 

Ai,„A,,„A„<,A„., ^ g^p |^_J_ j A,.h^(l + Jdec) 

This is consistency relation for the cosmological model with the slow-roll inflation. 
The meaningful cosmological should fulfill the condition 6 p^ 1. At present, the 
application of (|3.23|) is limited due to the unknown factor Areh- However, it could 



be possible to determine this value basing on the found value of m and the decay 
rate of the inflaton field. This issue requires however detailed studies, therefore 
we leave it to investigate elsewhere. 
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4. Modified inflationary spectrum and the CMB 

As it was shown in Introduction, the spectrum from the slow-roll inflation can 
be parametrized in the power-law form. However, the prior phase of a bounce 
should results with modification of this spectrum. The modifications of the pri- 



mordial scalar spectrum, were investigated in the numerous papers p2|, ^, P4 
However, the studies were performed only when the evolution of the scalar modes 
hold the classical form. Within loop quantum cosmology, not only the dynam- 
ics of the background is modified but also the perturbations |^ ^. In case 
of the tensor modes (gravitational waves), the form of these modifications was 
studied in details ^7\, ^ ^, ^. Based on this, the spectrum from the shark fin 
scenario considered has been recently found in p. In this paper both quantum 
corrections to the background as well to the perturbation part were taken into 
account. The case of the scalar modes is however more problematic. It is because 
of the issue of quantum anomalies of the algebra of constrains. In case of the 



so-called inverse- volume corrections, this problem has been resolved |^. How- 
ever, in case of the holonomy corrections the anomaly free equations are still not 
available. Some preliminary attempts to investigate LQC effects on the scalar 



power spectrum were performed in Ref. [B2|. However, the evolution of modes 



was treated classically and the quantum effects were introduced by the influence 
on the matter part. Also so attempts to derive holonomy corrected equations on 



scalar modes were performed in Ref. |^. However, the authors neglected the 



issue of the anomaly free algebra of constraint. Therefore derived equations can 
drive the system out of the surface of constraint and lead to erroneous predic- 
tions. Therefore, the systematic analysis of the anomaly freedom in case of the 
holonomy corrected scalar perturbations remains to be done. 

Before the details of modifications of the scalar power spectra will be avail- 
able, it is advisable to perform the phenomeno logical analysis of the possible 
impact of these effects on the CMB spectrum. This issue of impact of the LQC 



effects on the CMB spectrum was also discussed in Ref. [^ 0, |3^. However, 
here we perform quantitative analysis in contrast of the much more qualitative 
discussion in the cited papers. We assume that the scalar power spectrum takes 
a form 

Vs{k) = A{k,k,)Aj-j . (4.1) 

This is, in fact, the standard inflationary spectrum modified by the additional 
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prefactor A{k, A;*). The bounce-factor A(A;, A;*) can be written as 



A{k,K 



sm 



3k 

2k, 



3k 
2fc, 



(4.2) 



which is the simphfied form of the expression found in @]. The fc* is a parameter 
of the model and its interpretation will be discussed later. The factor A{k,k^) 
reflects typical modifications which appear in the bouncing cosmology. In the 
UV limit, limfc^oo A(fc, fc^,) = 1, therefore the spectrum ( [I.l|) is recovered. In 
turn, in the IR limit, \imk^oA{k, k^) = 0, and the spectrum is suppressed. This 
behavior of the power spectrum is typical for the bouncing cosmologies. The two 
effects of the bounce are transparent: suppression on the low k and the additional 
oscillations. In Fig. ^ we show function A defined by equation ( [4.2| ). Instead 
of using the wavenumber k we have translated it to the corresponding length 
A = ^, respectively A* = |p. 

In Fig. ^ we also show the 
function 

A(A,A.)^l + ^£. (4.3) 

This function measures the mod- 
ification due to the oscillations 
for A/ A* < 1. At A/ A* ^ 1 the 
spectrum becomes suppressed. 
In the bouncing cosmology the 
length scale A* can be related 
with the scale of horizon at the 
beginning of inflation. This is- 
sue was discussed in details in 
. Therefore if the present value 

of the scale factor is equal oq = 1 (as used in this paper), we have k^, ~ atHi 
where Oj is the value of the scale factor at the beginning of inflation and Hi is 
the value of the Hubble factor at the same time. Therefore if k^, and H^ could be 
measured, the total increase of the scale factor, from the beginning of inflation 
till present, can be determined. The value of fc* and respectively A* which is a 
scale of suppression in the spectrum is the crucial observational parameter of the 
bounce. In this paper we make an attempt of determining this value based on 
the observations of the CMB. 




0.001 



Figure 4: Plot of function A defined by equation 
(|4.2[) (solid line). The dashed line represents ap- 
proximation 
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As mentioned earlier, beside the effect of suppression, also oscillations of the 
spectrum are predicted within the bouncing cosmologies. This effects is much 
weaker that suppression, however is present also on the much smaller scales. 
This is important from the observational point of view. Namely, the length scale 
A* = 1^ can be much larger than the present size of horizon {k/k^, ^1). Then, 
the effect of suppression would be inaccessible observationally. However, some 
oscillations are still present on the sub-horizontal scales. Of course the amplitude 
of these oscillations decreases while k/k^ ^ 1. If the scale A* is however not 
much higher than the size of horizon, the effect of sub-horizontal oscillations 
could be quite significant. The oscillations in the primordial power spectrum 
translate into the additional oscillation in the spectrum of the CMB anisotropies 



(see e.g. [0). For the small multipoles, this subtle effect can be dominated by 
the contribution from the cosmic variance. However, for the larger multipoles this 
effect can dominate. At these scales, improvement of the instrumental resolution 
are still possible, what gives the chance to, at least, put a stronger constrain on 
these effects. 

In this section we confront the spectrum (^TTI) with observed anisotropies of 
the cosmic microwave background radiation. We use the seven years of obser- 
vations made by the WMAP satellite [jl|. In the numerical calculations we use 
the pubhcly available GAME code and Markov Chain Monte Carlo (MCMC) 



package CosmoMC |Q together with the CosmoClust code [^ for computing 
the Bayesian evidence. The codes were suitably modified to investigate the spec- 
trum ( [4.11 ). In computations, we take the standard cosmological parameters as 
follows 

{Ho, rtbh"^, n^h^, t) = (70, 0.0226, 0.112, 0.09) (4.4) 

and the pivot scale ko = 0.05 Mpc~^. 

In Fig. ^ we show spectrum of the CMB temperature anisotropies obtained 
based on power spectrum ([l.ip. The case A(fc, k^) = 1 (A; — )• 0) corresponds to 
the classical case with no contribution due to the bounce. The blue line corre- 
sponds to the best fit case. In this case, the modulations on the low multipoles 
are well reproduced. This is due to the oscillations in the primordial power spec- 
trum ( [4.1|) . This suggests that the effects of oscillations in the primordial power 
spectrum can be indeed studied basing on the CMB data. Perhaps the anoma- 
lous behavior of the CMB spectrum at / ^ 20 and / ~ 40 could be also explained 
by the oscillations within the bouncing scenario. However, not basing on the 
parametrization employed in this paper. The amplitude of oscillations on the 
lower scales must be higher than predicted by our model. 
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Figure 5: Spectrum of the CMB anisotropy. 

We also find confidence intervals for the parameters of the model, namely 
on As, Us and k^,. In these computations we take into account the temperature 
anisotropy data (TT spectrum) as well as the polarization data (TE and EE 
spectra). We neglect a contribution from the tensor modes putting Vt = 0. We 
show the obtained confidence intervals in Fig. ^ As it can be seen from Fig. ^, 
the parameters Ag are Ug are constrained from the both sides. Based on the fit 
to the WMAP data we find 



n^ 



0.97 ±0.07, 



As = 2.1- 10"^ ±0.1- 10"^ 



These results are in agreement with ( |3.2D and p.l| ). However it must be pointed 
out that they were computed at the different pivot scales. 

The parameter A;* has the upper constraint, however it is unbounded from 
below (large scales). It could be expected, since there is no observational data 
on the largest (greater than the Hubble radius) scales to get the upper constrain 
the parameter. Nevertheless some particular value of k^ is privileged what leads 
to the peak in the probability distribution. Based on the fit to the WMAP data 
we have obtained following values of this parameter: 

k^ = 1.7 ■ 10"^ ± 0.8 ■ 10"^ [Mpc-^]. 
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Figure 6: Constraints for the parameters As, Ug and k^. 2D plots: solid lines show 
the 68% and 95% confidence intervals. ID plots: dotted lines are mean likelihoods of 
samples, solid lines are marginalized probabilities. 



The length scale corresponding to fc^, is equal to 

27r 
A, = — ^ 4 ■ 10^ Mpc. 



(4.5) 



5. Big Bang vs. Big Bounce 

In this section we compare the model with suppression with the standard in- 
flationary model. The suppression appears generically within the Big Bounce 
cosmology. In turn, the slow-roll inflation in the standard Big Bang scenario 
does not lead to any suppression. While the suppression introduces a new length 
scale, the model with suppression has one more parameter in comparison with 
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the standard case. The considered models are: 

Hi - The slow-roll inflation within the Big Bang cosmology (spectrum given 
by Eq. |1.1| ). This model has two parameters A^ and rig. 

H2 - The slow-roll inflation within the Big Bounce cosmology (spectrum given 
by Eq. [4.1|). This model has three parameters A^, Ug and A;*. 



In the Bayesian approach to model comparison the best model has the largest 
value of the so-called posterior probability in the light of data, which is defined 
in the following way PD| : 

PiH,\D) = PJ^WPW (5.1) 

The Hi stands for considered model and D denotes data used in analysis. P{Hj) 
is the prior probability for the model under investigation, which should reflect 
all information which we have about it before the analysis with the data D, that 
comes from theoretical investigations, or from analysis with other data sets. In 
particular, if we have no foundation to favor of one model over another one, 
which is usually the case, we take equal values of P{Hi) for all considered mod- 
els. P{D\Hi) is the marginalized likelihood function over the allowed parameters 
range, which we called evidence and is given by 

Ei = P{D\Hi) = j deL{e)P{§\H,). (5.2) 

The 6 denotes vector of model parameters, L{6) is the likelihood function for 
considered model and P{9\H.j) is the prior probability distribution function for 
model parameters. 

It is convenient to consider the ratio of models probabilities, which is reduced 
to the evidence ratio (so called Bayes factor) when all considered models have 
equal prior probabilities: 

B,, = |i. (5.3) 

Their values give us information about the strength of evidence in favor of better 



model [^: if < In 5 < 1 we could not give conclusive answer, if 1 < Ini? < 2.5 
there is weak evidence, if 2.5 < Ini? < 5 the evidence is moderate, and for 
In i? > 5 evidence is strong. 

The values of evidence for two alternative models of primordial perturbation 
spectrum was calculated with the help of CosmoClust code, which was introduced 
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by [|^ as a part of CosmoMC code. We have based on anisotropy (TT) and 
polarization (TE, EE) data from the WMAP satelhte. In computations, we have 
neglected the contribution from the tensor power spectrum (we set Vt{k) = 0). 

We assume that models are equally probable (P(ifi) = -P(-f^2) = 1/2)- We 
consider flat prior probability distribution functions for unknown parameters in 
the following ranges: A, e [1.5 ■ 10~^ 5.5 ■ 10~9], n, e [0.5, 1.5], K G [10-^ lO"''^]. 
The value of logarithm of the Bayes factor which was obtained in the analysis, 
i.e. 

ln{Ei/E2) = In B^ = 0.2 ± 0.6, (5.4) 

does not give a conclusive answer. The data was not informative enough to 
distinguish these models. Therefore, in the light of the recent WMAP data the 
Big Bang and Big Bounce cosmologies are indistinguishable. The Big Bounce 
predictions are not in conflict with the observational data. Moreover, beside the 
fact that the Big Bounce model has one more parameter k^,, the obtained evidence 
is comparable with the Big Bang case. 

The above result was obtained with use of the CosmoClust code which bases 
on the nested sampling method ^^. This method was applied also in the Cos- 



moNest code [^. The computations with use of CosmoNest gives lni?i2 = 
1.1 ± 0.2. Therefore a week evidence for Big Bang model is obtained. How- 
ever, the CosmoNest was designed only for the case of the unimodal likelihood 
functions. In turn, the CosmoClust code extends to the case of the multi-modal 
likelihood functions. As it is clear from the bottom right panel in Fig. |^, the con- 
sidered likelihood function (dotted line) is bimodal in the subspace k^,. The first 
peak is located at k^ ~ 1.5 ■ 10~^Mpc~^ while the second at k^ ~ 2.5 ■ 10~^Mpc~^. 
Therefore the results from CosmoClust are more relevant for our model. The 
CosmoNest samples only around the highest peak, neglecting the contribution 
from the smaller one. Because of this, the observed discrepancy between the 
CosmoClust and CosmoNest results appears. It is worth to note that, the similar 
model with suppression on the large scales was shown as an example of use of 
the CosmoClust code pD[ . The bimodality of the likelihood functions was also 



observed and applicabihty of the CosmoClust code to that cases was emphasized. 

The issue of constraining the bouncing cosmology with the observational data 

was raised before in literature. In particular, studies based on SNIa data, location 

of acoustic peaks in the CMB and constraints from primordial nucleosynthesis 



(BBN) were performed in Ref. [^4], ^. However, these cosmographic methods 



are inefficient in searching for the effects of the bounce. It is due to the fact that 
the factor — is extremely low at the energy scales covered with this method. Even 
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during the BBN, where Tbbn ~ 1 MeV, we have pbbn ~ 10^^°ppi. Therefore, 
if pc ~ ppi, we have — ~ 10~^° and the holonomy corrections in the Friedmann 
equation (|1.2|) are vanishingly small ^. Based on the method developed in the 
present paper, we reach pobs = "^ 0°''° ~ lO^^^ppi, what gives — ~ 10~^^ for 
p^ K. ppi. Therefore, sensitivity on the holonomy corrections was increased around 
10*^° times with respect to the BBN constraint. 

Based the the results presented in this section one can conclude that the Big 
Bounce is consistent with the observations up to energy scales ~ 10~^^ppi. In 
this region the Big Bounce and Big Bang cosmologies are indistinguishable in 
the light of the available observational data. The advantage of the Big Bounce 
model is however that the initial singularity problem is resolved and the initial 
conditions for the phase of inflation are naturally generated. 

6. Can we see the Big Bounce? 

The present value of scale A,, is crucial from the point of possible observational 
investigations of the Big Bounce cosmology. As it was discussed before, this scale 
overlaps with the size of the Hubble radius at the beginning of inflation. There- 
fore, it corresponds to the point of maximal displacement of the infiaton field, 
namely 0max- In this section we investigate how the variation of 0max influences 
on the present value of A*. Based on this, it will be possible to investigate the 
observational conditions on the bounce. 

In Fig. ^, the schematic illustration of the scalar field evolution near the 
place of the maximal displacement was shown. In this figure we have marked 
the discussed 0max value as well as the observed value 0obs = 2.9mpi. While 
= 0obs5 the modes of the present size Aq = 3.14 Gpc (pivot scale) were formed. 
Based on this, we can determinate what is the present size of the mode, which 
was equal to the Hubble radius at = 0max- The transition from = 0max to 
~ corresponding to the total amount of e-foldings from inflation, which can 
be decomposed as follows A^tot = ^N + A^obs- Here A'obs is the observed value 
which corresponds to the transition from (f) = 0obs to ~ 0. The number of 
e-foldings during the transition from 0max to 0obs can be expressed as follows 

^^=-^ V'^'^ = ^ ('^n.ax - 0L) • (6.1) 



•^The constraint from the BBN can be however more significant in case of the so-called inverse 
volume effects in LQC |43] 
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Based on this expression as well as on the Friedmann equation, the present 
value of A* can be expressed as follows 



A* — An 



mpi 



mpi 



''obs 



exp < 27r 



mpi 



27r 



"obs 



nip] 



(6.2) 



^bs 



where Aq = 3.14 Gpc and 0obs = 2.9mpi. In Fig. |] we plot function A^, 
given by ( |6.2| ). For comparison, we also show some relevant length scales. 

The first one is the Hub- 
ble radius Hq/c ~ 4 Gpc. The 
second is the distance to last 
scattering shell (LSS), -Dlss ~ 
14 Gpc. The last scale is the 
scale of suppression A* ~ 40 
Gpc obtained in Sec. |. If 
0max > 2.94mpi then the scale 
A^, is placed behind the scale 
of LSS. In such a case there 
is no chance to see the effect 
of suppression directly. It is 
because, the scale of suppres- 
sion is higher than the phys- 




time 



Figure 7: Schematic illustration of the scalar field 
evolution near the place of the maximal displace- 
ment. The ^max is a maximal displacement of the 

• It,- r 1, - 1 1 field. The d>ohs is the value of the scalar field that 

ical horizon of photons, released ^ 

1 . ,, 1 • ^- mi corresponds to the powers spectrum measured at the 

durmg the recombmation. i here- 

r 1 -r , r^ nA pivot scale Aq = 3.14 Gpc. 

tore only it 0max < 2.94mpi, 

there is a possibility to study the effects of suppression on the CMB. From the 
fit performed in Sec. ^ we got A* ~ 40 Gpc, what correspond to 0max ~ 2.97mpi. 
Based on this, the one particular evolutionary trajectory can be distinguished. 
However, one have to keep in mind that the probability distribution on the param- 
eter fc* was unbounded from below. Therefore the obtained value 0max ~ 2.97rripi 
could be seen rather as a lower constraint on 0max- As mentioned, in order to 
make the direct observations of the suppression possible, the value of 0max should 
be smaller than 2.94mpi. The observations suggest that this value is higher, what 
unfortunately exclude this possibility. Based on this one can however exclude 
some models, where the predicted value of 0max is not higher than 2.94r?2pi. This 
is in fact a case for the symmetric inflation as studied in Ref. [§|. The issue 
of constraining this model was preliminary discussed in 0. This is also still 
possible that the effect of oscillations can be observed. Perhaps it is even the 
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reason why the particular value of A;* was distinguished from the WMAP obser- 
vations. Namely, it was possible because the structure of modulations at the low 
multipoles was reconstructed, not because the scale of suppression was detected. 

We finish this section with ^ [Gpcl 
discussion of the observational 
constraint on the the parame- 
ter pc- In loop quantum cos- 
mology, total energy density is 
constrained by p < pc- At the 

stage of inflation where the present ^.^.. ?.".^.^lf.™^.i"ll'^/!^''l 

pivot scale structure were ini- 
tiated, the energy density is equal 
to 
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Figure 8: The present value of the scale A* as a 

function of (/>max- 



Based on this, we infer that pc > Pobs- Because pobs ^ Ppi, the observed constraint 
on the energy scale of the bounce is very weak. However, since pc ~ 1/7^, the 
constraint on the parameter 7 can be much stronger. Indeed, based on ( |1.3| ) we 
find 



7 < 1100. 



(6.3) 



The value obtained from consideration of black hole entropy 7 = 0.239 places well 
within the observational bound. The constraint ( p. 3D is quite strong, however it 
must be kept in mind that it is based on relation ( |1.3|) , which can be invalid. 



As discussed in Ref. |13[, the pc can be a free parameter, and then it would be 
impossible to put the constraint as ( |6.3|) . Therefore, more theoretical predictions 
regarding the phenomenological parameters as, pc, are still awaiting. 

To conclude, some models of the bouncing cosmology can be excluded based 
on the observations of CMB. It is based on the observational constraint on 0max- 
The direct observations of the bounce effects are however much harder to detect. 
As we have indicated, the effect of suppression cannot be used. It is because the 
scale of suppression was shown to be higher than the scale of horizon. The effect 
of oscillations gives a chance, however the effect is, in general, weaker and can be 
below the cosmic variance. It must be also pointed out that the discussed effects 
can be also predicted from the different models. Therefore the important task is 
to find the observable which enables to distinguish between the models. 
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7. Summary 



In this paper we have examined observations of the cosmic microwave background 
radiation as the potential probe of physics in the Planck epoch. We have based 
our considerations on the predictions of loop quantum cosmology. Within this ap- 
proach the initial singularity is replaced by a cosmic bounce. During the bounce, 
the universe reaches the maximal energy density pc < oo. In our studies, we 
have concentrated on the model with a massive scalar field. The advantage of 
this choice is that the phase of the standard slow-roll inflation is realized during 
the evolution. During this phase, the primordial perturbations are formed. Their 
properties can be investigated by observations of the CMB. Based on the recent 
results from the seven years of observations made by the WMAP satellite, we 
have determined some parameters of the model. In particular, we have found 
that inflaton mass m = (1.4 ± 0.5) ■ IQ-^mpi = (2.6 ± 0.6) ■ 10^=^ GeV. 

Subsequently, we have investigated the modifications of the primordial power 
spectrum due to presence of the bounce. The two main effects that were discussed 
are: suppression and oscillations of the inflationary spectrum. The suppression 
can very strongly modify the spectrum. However, the present scale of suppression 
was shown to be behind the Hubble radius. Therefore, the effect of suppression is 
not directly observationally available. Despite this, the present scale of suppres- 
sion A^, = 1^ ~ 40 Gpc was distinguished by the observations. This is because of 
the oscillations in the primordial power spectrum. The effect of oscillations is in 
general, more subtle and dominated by the cosmic variance. However, the present 
work indicates that the oscillations in the power spectrum can in fact explain the 
strong additional modulations in the spectrum of CMB on the low multipoles. In 
order to verify it, we have compared the Big Bounce model with the standard Big 
Bang scenario and showed that the present observational data is not informative 
enough to distinguish these models. In other words, the Big Bounce predictions 
were shown not to be in conflict with the observational data. Moreover, despite 
the Big Bounce model has one more parameter /c*, the evidence obtained in this 
case is comparable with the Big Bang case. 

We have shown that 0max is not lower than 2.97r?7,pi. Based on this, some 
models of the bounce, as the mentioned symmetric model, can be significantly 



constrained or even excluded. Assuming validity of (|1.3|) we gave an observational 
constraint on the Barbero-Immirzi parameter. We found that 7 < 1100. The 
corresponding constraint on pc was shown to be much weaker. 
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